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Objectives 


•  Advance  the  state-of-the-art  in  multi-resolution  algorithms  and  parallel  software. 

•  Develop  high-order  adaptive  discretizations  on  unstructured  meshes  for  the 
Navier-Stokes  equations  (mostly  incompressible)  at  high  Reynolds  number. 

•  Apply  algorithms  to  parallel  simulations  of  complex  phenomena  in  turbulence 
and  other  transport  phenomena. 


List  of  Accomplishments 

•  Development  of  a  new  hierarchical  specrtal  basis  constructed  appropriate  for  hp- 
finite  element  formulations  on  unstructured  grids  consisting  of  triangular  and  tetrahe¬ 
dral  subdomains.  This  basis  uses  non-symmetric  highorder  Jacovi  polynomials  and  is 
the  building  block  for  the  accurate  (spectral)  solution  of  partial  differential  equations 
in  complex  geometry  computational  domains. 

•  In  more  detail,  the  new  multi-resolution  basis  has  the  following  properties:  (1) 
Jacobi  polynomials  of  mixed  weights,  (2)  Semi-orthogonality,  which  leads  to  great 
sparsity,  (3)  Hierarchical  structure,  (4)  Generalized  tensor  products,  (5)  Mixed  order 
expansion  which  provides  great  flexibility  in  adaptive  discretizations,  and  (6)  auto¬ 
mated  Gaussian- Jacobi  quadratures  which  minimize  errors  in  complicated  geometries. 
The  scaling  of  the  eigenspectra  of  the  first  and  second  derivatives  has  been  studied 
in  detail  in  order  to  determine  the  conditioning  of  this  basis.  We  have  shown  that 
the  scaling  of  these  derivatives  is  identical  to  the  scaling  in  standard  tensorial  spec¬ 
tral  bases.  In  addition,  due  to  a  new  type  of  tensor  product  we  have  developed,  the 
computational  cost  in  function  evaluations,  squares,  derivatives,  etc.  is  similar  to  the 
tensorial  spectral  methods  used  for  structured  discretizations. 

•  The  accuracy  of  the  new  method  has  been  tested  both  in  two-  and  three- 
dimensions.  For  smooth  solutions  we  obtain  exponential  convergence  irrespective 
of  the  triangulization.  Based  on  this  new  formulation  a  flow  solver  (NEKTAR)  was 
developed  in  order  to  be  used  in  turbulence  simulations.  There  are  two  important 
new  features  of  this  code  that  make  it  superior  to  our  previous  spectral  element  flow 
solvers.  First,  it  is  fully  adaptive  as  h-  and  p-type  refinement  can  be  readily  ap¬ 
plied  using  the  hierarchical  structure  of  the  basis.  For  example,  each  edge  or  face 
or  the  interior  in  each  tetrahedron  could  have  different  modes  according  to  resolu¬ 
tion  demands.  Second,  it  is  more  efficient  as  more  effective  preconditioners  can  be 
constructed  using  multi-resolution  type  algorithms. 
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•  An  important  numerical  result  was  obtained  verifying  that  that  high-order  accu¬ 
racy  can  be  obtained  on  standard  unstructured  CFD  meshes.  Moreover,  unlike  low- 
order  discretizations  which  may  require  Delaunay  triangulizations  for  convergence, 
the  new  hierarchical  basis  can  be  used  in  conjunction  with  highly  distorted  meshes 
without  accuracy  degradation. 

•  Development  of  a  new  projection  formulation  of  hyperbolic  and  parabolic  prob¬ 
lems  using  a  discontinuous  Galerkin  formulation.  This  allows  L 2  expansions  to  be  used 
even  for  second-order  elliptic  problems,  which,  in  turn,  provides  a  great  discretization 
flexibility.  For  example,  appropriate  trial  bases  can  be  used  in  different  subdomains 
to  deal  with  singularities,  discontinuities  and  sharp  gradients  while  optimizing  their 
efficiency  locally  without  global  inter-dependence. 

•  Extensions  to  two-dimensional  hybrid  discretizations,  where  both  quadrilateral 
spectral  elements  and  triangular  elements  are  connected  together  to  provide  greater 
flexibility  in  discretization. 

•  Extensions  to  three-dimensional  hybrid  discretizations,  where  polymorphic  ele¬ 
ments  can  be  used,  i.e.  hexahedra  connected  to  tetrahedra  or  pyramids  or  prisms. 

•  A  detailed  study  of  the  spectral  convective  operator  on  unstructured  and  hybrid 
meshes,  with  bounds  and  semi-analytical  expressions  derived  for  the  growth  of  the 
maximum  eigenvalue  as  a  function  of  the  spectral  order.  This  is  important  in  choosing 
appropriate  time  stepping  algorithms  for  convection-dominated  phenomena. 

•  Development  of  a  Navier-Stokes  solver  on  hybrid  meshes  and  validation  verifying 
spectral  convergence/accuracy  for  prototype  solutions  on  complex  geometries. 

•  Development  of  an  Arbitrary  Lagrangian  Eulerian  (ALE)  formulation  for  spec¬ 
tral  methods  on  unstructured  and  hybrid  mehes. 

•  Development  of  a  new  computing  approch  in  simulation,  namely  steering  com¬ 
puting,  where  the  developed  hierarchical  spectral  methods  are  used  on  grids  changing 
continuously  in  time  by  following  an  adaptive  procedure.  We  have  developed  and 
performed  such  steering  computational  procedures  for  parallel  computers  in  conjunc¬ 
tion  with  interactive  graphics  in  order  to  introduce  a  new  mode  of  Direct  Numerical 
Simulation,  i.e.  the  dynamic  DNS. 


Personnel 

•  Faculty:  G.E.  Karniadakis,  Professor  of  Applied  Mathematics 

•  Post-Docs:  T.  Matushima,  S.J.  Sherwin  (part-time)  and  A.  Beskok  (part-time) 

•  PhD  Students:  I.  Lomtev,  T.C.  Warburton  and  C.  Evangelinos. 
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•  ICOSAHOM’95  •  IBM  •  MIT  •  Cornell  Theory  Center  •  Boston  University  • 
Virginia  Polytechic  Institute  •  Hong  Kong  University  •  International  Conference  on 
Parallel  Algorithms  and  Applications  (Wuhan,  China)  •  Institute  for  Scientific  Com¬ 
puting,  Chinese  Academy  of  Sciences  •  Institute  of  Mechanics,  Chinese  Academy  of 
Sciences  •  University  of  Tokyo  •  University  of  Notre  Dame  •  University  of  California 
Santa  Barbara  •  Wright  Patterson  Air  Force  Base  •  City  College/Levi ch  Institute  • 
AFOSR  Boiling  AFB  DC  •  Nuclear  Regulatory  Commission  •  National  Institute  of 
Standards  and  Technology  •  AIAA  Conference  New  Orleans  •  University  of  Rhode 
Island  •  University  of  Maryland/CTC  •  University  of  Cincinnati  •  Worchester  Poly¬ 
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•  First  AFOSR  Conference  on  DNS/LES  •  ASME  micro-Therm  Workshop  •  10th 
Domain  Decomposition  Conference. 

The  code  NEKTAR  has  been  disributed  to  more  than  one  dozen  Universities 
and  Laboratories.  Some  of  them  include  Cornell  University,  Penn  State  University, 
University  of  Wisconsin,  Imperial  College,  North  Carolina  University,  Florida  State 
University,  OAK  Ridge  Labs,  Nielsen,  Inc.,  Boeing  Inc.  etc. 
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